Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN KTgX style file v2.2) 



Finite source effect on the polarization degree induced by a single 
microlens 



Hiroshi Yoshida 1 * 

in 

f**"^ , 1 Department of Physics, Fukushima Medical University, Fukushima-Citv 960-1295, Japan 

o 

£jjrj. 5 February 2008 



1 INTRODUCTION 



ABSTRACT 

We investigate the effect of a single microlens on Stokes parameters. Semi-analytical formulae 
of the microlensed Stokes parameters are derived. The formulae not only reduce the double 
integrals in the estimations of those quantities but can also be approximated to a useful form 

■ in the bypass case. By using our formulation, we show that a combination of polarimetric 
, data with photometric data enables us to estimate not only the finite source effect but also the 

^"T 1 ' direction of the microlens motion. 
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$_l ' produced by a compact object m ttie toreground, such as a star or a planet. tFaczvnsKil !19Kcl) sug gested mat microlensing events could be 
r/3 , detected by a succe ssive observation of stars in LMC, SMC or our Galactic bulge, and since 1993 iAlcock et alll993llAubourg et alll993l: 
C$ ■ ludalski et al]|l993h many observational groups have reported a substantial number of microlensing events. According to IPaczvrrskjT the 
observations of microlensing events provide information on the mass distribution of invisible objects in galaxies; further, they impose some 
restrictions on the dark matter. 

The observation of a microlensing event provides a light curve of the source, which is a fundamental observable in this event. From 

■ analysis of the light curve, it is important to estimate the microlens parameters such as the mass M of the microlens object, the distance 
parameter £(= Dol/ Dos, where Z?ol,-Dos are the distances from the microlens and from the source to the observer, respectively) and 
relative velocity V of the microlens orthogonal to the line of sight. However, it is usually very difficult to obtain the value of each parameter 
because these parameters are coupled with each other in observed quantities, e.g. the event duration t = 2tv = 2r-E/\V\ and the normalized 
minimum impact parameter uo = b/rv, where te = (4GA'/DosC(l — C) /c 2 ) 1//2 is the Einstein ring radius on the lens plane. 

In order to decouple the degeneracy in the mic rolens parame ters, several authors have proposed useful effects associated with microlens 
event s , namely, the para llax effect on the light cu rve lGouldll992l) . and the accurate astrometry of the lensed source hioie. Novikov & Polnarevl 
Il99.4 iMivamoto & Yoshiilll99.4 IWalkejTl995h . The former effect enables us to measure the Einstein ring radius projected on the ob- 



The microlens effect is a phenomenon in which the flux emitted from a so urce (star) varies with time due to the gravitational lens effect 
produced by a compact object in the foreground, such as a star or a planet. |Paczvhskilfl98fitl suggested that microlensing events could be 
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2003; 



lMaolll999tlAlcock et alJhood bond et alj 



teKI — C))- This effect has been observed in a large number of 



200 ll; iMao et allhool Ismith et aljEool ISmith. Mao & Paczvnskil 



l2003bh . The latter effect provides the angular Einstein ring radius 



0e(= t^/Dol) through the observation of the fine motion of the lensed source. In order to detect this effect, we need the /^as level accuracy 
in the astrometry. Hence, some satellite projects are being planned for carrying out such observations. These two effects are expected to 
impose restrictions on the determination of the microlens parameters. In particular, the detection of both the effects in the same event enables 
us to decouple the degeneracy in the parameters M, £ and \V\ as M — c 2 Ye#e / '4G ', C = Fe/(-Dos$e + Pe) and |V| = C^osfe/te, 
respectively. 

To date, certain events that involved the estimation of the microlens mass without the high-precision astrometry have b een reported 
ISmith et aj2003bllYoo et all2004lpark et alfeoollGhosh et aj2004l ; |jiang et alfcooA . In four of these events, the authors ISmith et all 
IYoo et all IPark et ali Lliang et alii not only detected the parallax effect but also estimated t he angular Einste in radius by combining an 
empirical relation between the angular size 9* and the colour-magnitude of the source given bv lvan Bellel Il999h with the finite source effect 
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p* = Ot/O-E on the light curve Iwitt & Maolll994l : iNemiroff & Wickramasinghelll994[) . The finite source effect is relatively easy to be 
detected in the transit case, where the lens transverses the face of the source projected on the lens plane; this is because the light curve of 
an extended source is significantly differ ent from that of a point source when the normalized impact parameter u is such that u < p, (e.g. 
Alcock ct al ll997tlSoszvnski et alhOQlh . The above events that are successful in decoupling the degeneracy in the parameters are examples 
in the case. On the other hand, in the bypass case, where the lens never transverses the fa ce of the sour c e, it is difficult to detec t the effects 
because we can fit the light curve using appropriate parameters in a point source model iOoul dll994h . Idould & Welch] Jl995» found that 
in the bypass case the finite source effect is not detectable by a single-band observation but can be detected by a two-band observation, 
especially in t he optical and infrared bands. 

Recentlv. lsimmons. Willis & Newsam|ll995ah and lsimmons. Newsam & WilliJh 995th pointed out that t he finite source effect also ap- 
pears in another aspect of the microlens effect. They found that the Chandrasekhar effect IChandrasekhaJl96Ch is enhanced by the microlens 
effect and that this enhancement causes the polarization degree of the source to vary with the motion of the microlens. They also found that the 
polarization degree has a double peak in the transit case, while it has a single peak in the bypass case. They presented a numerical calculation 
of the Stokes parameters with the microlens effect; however, they did not present the analytic al formula, although the size dependence of the 
polarization degree helps us to decouple the degeneracy in the microlens parameters. Before lsimmons et alllSchneider & Wagoned Il987h 
also discussed the gravitational lens effect on the polarization degree as a tool for identifying a lensed supernova. In their paper, by using a 
simple approximation of the magnification factor 1 for a point lens, they presented a semi-analytical formula of the polarization degree and 
the magnification factor to demonstrate that the size expansion during the supernova explosion yields variations in these quantities. However, 
they did not discuss the benefits from the polarimetry of the microlensed source at all. 

In this paper, we focus on the finite source effect on the Stokes parameters and the polarization degree of the microlensed star. Further, we 
will also discuss the possibility that the polarimetric observation of the lensed source decouples the degeneracy in the microlens parameters. 
Throughout this paper, we consider a single microlens. The outline of this paper is as follows. In Section 2, we provide a brief review 
of the microlens effect. In Section 3, we present semi-analytical formulae of the Stokes parameters to discuss the finite source effect on 
the polarization degree. Then, we provide an approximation of the formula with regard to the bypass case. In Section 4, we consider the 
extraction of information concerning the finite source effect from the polarimetric observation of the microlensed source. Finally, in Section 
5, we discuss some of the possibilities of detecting the polarization induced by the microlens effect. 



2 MICROLENS EFFECT ON THE LIGHT CURVE 

When light rays from a distant source pass near a compact object (microlens) with mass M, they are deflected by an angle AGM/<?b, where 
b is an impact parameter. Let Dos, Dol be the distances from the source and from the microlens to the observer, respectively. Let 9s be 
the angular position of the source that would be observed without the microlens. In this situation, the lensed light (hereafter referred to as 
'image') is observed at the angular position 6, which is obtained as a solution of the following lens equation: 

s = e-$o, (i) 



where 9e = y/AGM(l — C)/C c2 Dos is the angular Einstein ring radius. The distance parameter £ is given by Dol / Dos ■ In equation 0, 
let us introduce the dimensionless variables u = 9s /9b and x = 9/9e and rewrite the lens equation as follows: 
x 

u = x--—. (2) 

\x\ 2 

In general, the equation has two solutions, x±, and they give the image positions of the source. One of the images is gravitationally magnified 
by factor A+, while the other is magnified by factor A- (demagnified): 



2 \uVu 2 + 4 

However, we cannot resolve each image because the typical separation angle between the images is ~ 9e < 0.1 — l(M/Mo) 1//2 mas. 
Hence, we observe them as a single image with the following magnification factor: 

A p (u)=A + (u)+A^(u) = ^±T7- (4) 

Uy/u 2 + 4 

Because the microlens moves with a velocity V relative to the line of sight, the normalized distance u from the microlens to the source 
varies with time as u = {ul + r 2 } 1/2 , where uo is the normalized minimum impact parameter of the microlens, r = (t — to) /is is the 
time normalized by the time-scale of the microlensing event (fE = Dol#e/|V|) and to is the time at which the magnification factor has a 
maximum value. Consequently, the microlens motion yields the variation of the magnification factor. 



According to lSchneider. Ehlers & Falcclll992h . we use the term 'magnification factor' instead of 'amplification factor' 
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source 




Microlens 

Figure 1. Geometry of microlens and source: u and w are the vectors to the centre and to a point marked as '+' at (p, <j>) of the source projected on the 
lens plane from the microlens, respectively. i9 and a are the angles between V and u and between V and the x-axis defined by an observer, respectively. 
Therefore, ip b s = + $ + a is the angle between the x-axis and the line joining the marked point and the centre, no is the minimum impact parameter. 



By taking into account the finite size of the source, some authors pointe d out a remarkable difference between the light curves of the 
point source (the standard microlens) model and of the extended source model Iwitt & Maoil994LlNemiroff & Wickramasinghell994^ . This 
effect is termed the finite source effect. In general, the lensed flux is given by 



S/(p.|u) 



O'ipdp 



7(6> E p, <t>)A p (w(u, p, 4>))&<j> = A(p*\u)Sio, £ro 



9'ipdp 



(5) 



where I(r, <f>) is the surface brightness and A(p*\u) denotes the total magnification factor. In the above equations, p denotes the normalized 
radial distance to a point on the lens plane from the centre of the projected source, while (j> denotes the angle between u and the line joining 
the point and the centre (Fig.0. Furthermore, p* = #*/#e (#*: the angular radius of the source) is the normalized radius of the source, and 
w is the normalized distance from the microlens to the point, given by w(u, p, ( 

I(8ep), is given by 



y^U 2 + 2up COS (f> + 0" 

The total lensed flux from the extended source with axial symmetry in the surface brightness, I(r, <j>) 



Si(p t \u)=40£ 



pdpl(6vp) 



(u - /9) 2 n(7r/2; k, n) + 2F(tt/2; k) 



(6) 



{u + p)y/( U -p) 2 +4 

where F(ty/2; k) and II(7r/2; k, n) are the complete elliptical integrals of the first and of the third kinds, respectively iGradshtevn & Rvzhild 



see also Appendix A). Further, 

4wp 



A: 



4/i 



(7) 



(u + p) 2 ' " Y(^-p) 2 + 4' 

were introduced by Iwitt & Mad where the total magnification factor of the flux from the source with a constant surface brightness was 
estimated. An equivalent expression of equation J5J was also presented bv lwidfl99l , which enables us to numerically calculate the total 
fluxes in various bands. Using this expression, we can obtain the total magnification for the source with the limb-darkening surface brightness 



I(0ep) = Io 



as follows: 



A(p.\u) = 



1-Cl 1-4 1- 



p- 



7rp 2 (l- Cl /3) 



dpp 



1 - ci 1 - 



1- P - 

Pl 



(it - p) 2 n(7r/2; k, n) + 2F(tt/2; k) 

{u + p)yJ{u-p) 2 +4 



(8) 



(9) 



where ci is a parameter that depends on the wavelength of the light. 

In Fig. [2] we present the contours of A{p*\u) on the it/p*-p» plane in panel (a) and the u/p* dependence of A(p*|u) in panel (b). 
Both panels clearly show that the behaviour of A(p,\u) inside the source is very different from that outside the source. Outside the source, 
the contours (Fig.|5Ja), particularly, for low p*, are almost straight lines with a slope of p*/u, which shows that the magnification factors 
depend only on u. On the other hand, the contours inside the source are rather flat, which means that the factors inside the source have little 
dependence on u. 
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Figure 2. Magnification factor: contours of the magnification factor (a) and the u/p, dependency of the magnification factor (b). In our estimation, the linear 
limb-darkening model I(0ep) = ^o{l — ci[l — (1 — p 2 / P*) 1 / 2 ]} is adopted, where ci = 0.64. 



3 MICROLENS EFFECT ON THE STOKES PARAMETERS 
3.1 Stokes parameters 

Ichandrasekhail J 196(1 numerically calcul ated the effect o f electron scattering in the atmosphere of a star on the light emitted from the star 
and found that it leads to polarization (the lchandrasekhail effect). According to lSchneider & Wagoned fl987i) . we can assume that the Stokes 
parameters Q, U, V at each point (p, (j>) on the surface of the star with axial symmetry are expressed as follows: 



Q(<M = I C2 ^1 - ^ 1 - ^ j , U = V = 0. (10) 

Normally, we cannot resolve the light from each point on a star except for the Sun. We only observe the total flux that is obtained by 
integrating the lights over the surface of the star. In this case, the total Stokes parameters, Sq, Su and Sv, are also defined as those integrated 
over the surface. In integrating the Stokes parameters over the surface, it must be done based on the same reference frame. In this regard, we 
must fix x- and y- axes such as in the West and North directions, respectively. As seen in Fig.Q we define the angle between a>axis defined 
by the observer and the line joining the centre of the star and the point (p, <j>) as <p h B - We thus define the total Stokes parameters as follows: 

Sq 

Su | = / #Epd/5<2(6> E p) / difi ohB | -sin2i£obs | • (11) 
Sv 

In this case, we can easily deduce that Sq = Su = Sv = because of the axial symmetry. Therefore, the light from th e usual circular star 
is unpolarized. Specifically, we cannot detect the Chandrasekhar effect in a star without any symmetry breaking. In fact, iKempetalj] 19831) 
detected this effect during an eclipse in a binary system. 




3.2 Microlensing on polarization 

Isimmons et alj ll995ah pointed out that due to the symmetry breaking by the differential gravitational magnification, the cancellation of 
polarization between symmetrical points on the surface is incomplete, and the flux is totally polarized as a result. Moreover, the magnification 
of light from each part of the surface also varies with the microlens motion. Consequently, the polarization varies with time. lsimmons et all 
ll995bl) also numerically calculated the polarization degree and found that the time profile of the polarization has a single peak in the bypass 
case, while it has double peaks in the transit case. Unfortunately, it is difficult to understand the finite source effect on polarization, especially 
the size dependen ce of polarization, only on the basis of the numerical calculation. 
According to lsimmons et alj Il995bh . the lensed Stokes parameters are given by 

Sq ^ - e 2 E pdpQ(8 E p) r d^obs ( cos2 ^° bs )A p (w(u,p, <j>))- (12) 
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Figure 3. Polarization degree: Contours of the polarization (a) and the it/P* dependency of the polarization degree (b). In this estimation, the limb-darkening 
model Q(#ep) = 7q c 2[1 — (1 — P 2 /p*) 1//2 ] is adopted, where C2 = 0.04. In panel (a), the contour spacing is 0.05 percent. 



The relation between <p bs and <f> is (fobs = (p + iS + ce, where # = — tan -1 (uo /r) and a is the angle between the direction of the microlens 
motion V and the a>axis (see Fig.0. From this relation, we can integrate the inner part of equation i 12t with respect to cj> and express the 
Stokes parameters as follows: 



Sq \ n , n / cos2(tf + a) 



(13) 



D(p t \u) 



2 J dp | 



i-4 

pi 



(u-pf 



2 p(u + p)yj{u - p) 2 + 4 



n(7r/2;fc,n) 



8m 2 p(m + p)y / (M-p) 2 + 4 



8u 2 p 



(14) 



where E(tt/2; k) is the complete elliptical integral of the second kind; k and n are given by equation 0. Hence, we obtain the polarization 
degree as follows: 



V S Q + S 'u _ D(p*\u) 



(15) 



Si(pt\u) Si(p*\u) 

Using equations ^6} and i 141 . we can reduce the estimation of the polarization degree from a double integral to a single integral in equations 
J5j and I12i . In the numerical estimation, we can use some subroutines (e.g. GNU Scientific Library 2 ) for the complete elliptical integrals in 
the these equations. 

Fig-0shows the contours of polarization degree on the u/p*-p* plane in panel (a) and the u/p* -dependence on the polarization degree 
in panel (b). From panel (a), it is clear that for p* < 0.05, the contours are parallel to the p*-axis at least up to uj p, < 10. This means 
that the polarization d egree depends only on tt/p* for p* < 0.05 and u/p* < 5, which is also shown in panel (b). As mentioned in 
ISchneider & WagoneJ. we also confirm that the polarization degree has the maximum value P m ax ~ 0.875 per cent at u/p* ~ 0.96 for low 
p*. On the other hand, for large p„, it is apparent from panel (b) that the polarization degree strongly depends on the s tar's radius. From panel 
(b), furthermore, we observe that while for u/p, < 1, the polarization is almost proportional to (it/p*) 2 , as shown bv lSchneider & Wagoned 
whereas for uj p* > 1, the behaviour of the polarization is not as simple as argued by them. This difference arises due to the fact that 
ISchneider & Wagoned adopted an approximation to the magnification factor for the point lens, while we refrained from doing this in this 
stage. Our expression (eq. 1141 ) is the exact form for the single microlens. 



2 http : / /www .gnu.org/software/gsl/index. html 
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Figur e 4. A Comparis on of approximations to the polarization degree: The left panel shows contours of the relative errors of an approximation given by 
ISchneid er & Wagonej] to the polarization degree estimated by using equations |S| and 1141 numerically. The right panel shows contours of the relative errors 
of our approximation i 1 81 to the numerical estimation. The contour spacing is 1 per cent. 



3.3 Polarization degree in the bypass case 



The use of equations and 1 1 3i . besides reducing the double integral in equations ^5) and 1121 . is also advantageous in obtaining an 
approximation of the polarization degree in the bypass case . Hence, in this sub section, we focus on the polarization degree in the bypass 
case. By using equations (A1)-(A3) given in the Appendix o f lwitt&MaoNl994 . we can expand the complete elliptic integrals with respect 



to k and n in the bypass case (p, < Uo and p» < 1) 

60 - 28d 



We can thus approximate formulae ^6} and i 141 as follows: 



Si(p*\u) = A p (u)Sio 



1 + 



u 2 + 1 



15 



5ci (u 2 + 2) (u 2 + 4) 2 



D{p*\u) 



7c 2 A p (u)Sio 
5(l-ci/3) ( u 2 + 4) 2 



1 + 



304 (it 6 + 3m 4 + 6u 2 + 5) 
147 (u 2 +2) (u 2 +4) 2 



(16) 



(17) 



While the finite source effect p* is a secondary effect in the magnification factor, it is the primary effect in the polarization degree. Thus, we 
observe that the polarization evidently originates from the finite source effect. 

Using the above equations, we can express the polarization degree in the bypass case as follows: 



P(p.|u) 



7c 2 



(p.A 



5(l-ci/3) (u 2 +Af 



1 + 



304 m 6 + 3u 4 + 6u 2 + 5 60 - 28ci 



u 2 + 1 



(18) 



147 (u 2 + 2)(u 2 + 4) 2 15-5ci (u 2 + 2)(u 2 + 4) 2 

Equation J 1 XI is similar to the one derived bv ISchneider & WagoneJ : however, our approximation seems more useful even in the region 
u > 1. Fig.[4]present a comparison of approximation and theirs. Although both approximations deviate from the full expression FBI with a 
factor 1.5 — 2 at w/p* ~ 1, the relative error of our approximation is l ess than 10 per cent fo r 1.2 < w/p* < 10 and p» < 0.1, whereas the 
region with the relative error less than 10 per cent of the one given bv lSchneider & Wagonej exists for 1.6 < u/ p, < 10 and p„ < 0.05. 
In the bypass case, the polarization degree has a maximum value P max at r = 0, as follows: 



7c 2 



80(1 - ci/3) Vmq 



1 + 



871 



+ 



5880 20(1 - ci/3) 



(-) : 



(19) 



for p, ^ Mo < 1 and 1.2 < wo/p* < 5. For ci = 0.64 and C2 = 0.04 ^Schneider & Wagonej) . we obtain the maximum value of the 
polarization degree, P max ~ 0.12(2p<, /u p) 2 per cent. This value is very small, but not always undetectable. In fact, a polarization degree 
of this order was detected bv lKemp et all who obtained the polarization degree ~ 0.01 per cent during an eclipse in a binary system. 



3 In the microlensing events reported so far that exhibit the finite source effect (transit cases), the assumption p* < 1 is a good approximation lAlcock et alJ 
1997). 
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4 EXTRACTION OF THE FINITE SOURCE EFFECT FROM POLARIZATION 

In the previous section, we investigated the polarization by a microlensing event. As mentioned in Introduction, combining the finite source 
effect with the parallax effect enables us to decouple the degeneracy in the microlens parameters. In this section, we shall consider the method 
of extracting information on the finite source effect from the polarization. 

4.1 Variability of the polarization degree 

Isimmons etai]ll995bh showed that the polarization degree varies with time by a microlensing event and that the time profiles (variabilities) 
are classified into two cases, i.e. transit case and bypass case. We thus discuss the finite source effect on the polarization degree separately 
for each case. 



4.1.1 Transit case 



In the transit case, the polarization de gree has a double-peak profi le. Therefore, we can measure the peak times £1,^2 (> ti) from the 
polarimetric observation. According to lSchneider & Wagoned ll987t) . it is theoretically known that the polarization degree has peaks near 
it = 0.96p* . Thus, by combining At = £2 — ti with the analysis from the light curve (uo, te), we can estimate p* as follows: 



1.04 



*H0 



(20) 



4.1.2 Bypass case 

Unfortunately, the polarization has only one peak in the bypass case. Hence, we cannot measure At, which is possible in the transit case. 
However, we may be able to use equation 1 1 81 to estimate p,. Provided that p» has a small value and 1 < u/p* < 5, we can approximate 
the equation as follows: 



P = 



U \ ti / 



^c2pl 



80(l-ci/3) : 



871 



+ 



5880 20(1 -ci/3) 



Using the same approximation, we can deduce A(p,\u)/A p (u) from equation 1 1 6t as 



a{pAu) 

A p (u) 



= 1 + 



7_ 

40 



20(1 - ci/3) 



P* \ * -, c 

— =l+-r, C = 



J_ 

40 



20(1 - ci/3) 



(21) 



(22) 



Therefore, we can obtain the values of a, b, c from the data fitting of the light curve and the variability of the polarization degree in the same 
event. Consequently, the coefficients a, b, c yield p„, Ci and c 2 as follows: 



7600a 



294,, , 15(1476 - 233c) 

(0 + C), Cl = - and C2 = — ; r 

95 v ; 10296 -871c 7(10296 - 871c) 



(23) 



4.2 u/p* dependence of polarization for small p* , u 

As mentioned in subsection 3.2, for small p* and u, the polarization degree depends only on u/p,. From equations ^6) and 1 1 -41 . as well, 
we can prove the validity of this statement. In these equations, the arguments of the complete elliptical integrals, k and n are approximated 
as k 2 w n = 4t/(l + f) 2 for p, u < 1 (but < t = p/u < 00). By using the relation (u - p) 2 n(7r/2; k, n) ~ (p + u) 2 E{tt /2;k) 
derived from equation <A10> in Appendix A, we find that F(ir/2; k) terms are the leading terms in the integrands of equations {6} and < 141 . 
Therefore, we can approximate these equations as follows: 



&(p*|n) 
and 

D(p,\u) 



4S7 



tt(1 - ci/3) pi 

4C2<SlO u 

tt(1-ci/3) pi 



dt i 1 - ci 



1-y/T- (ut/ P *y 



tF(n/2;k(t)) 
1 + t ' 



(If 



1 - V 1 - (^/P.) 2 



(3f 4 + 2f 2 + 3)F(tt/2; fc(t)) 
87(1+7) ' 



(24) 



(25) 



From these equations, we have proved that P(p*]u) — D(p*\u) / S 1 (p,\u) depends only on u/p, for p» , u <C 1. 
Hence, we define the normalized Stokes parameters Sq and Su as 

Sq = Sq/Si = P(p* \u) cos 2(i? + q) and St/ = Sv/Si = -P(p*|«)sin2(i? + a), (26) 

respectively, where •& = — tan -1 Uo/t and a is the direction of V to the a>axis. Then, we can draw QU plots for some values of uo/p*. 
We have presented various QU plots for various uo/p* values in Fig. [5] where for all the curves, we employ the same values of p* and 
q; p* = 0.001 and a = 0. These curves are very different from each other. In particular, the curves in the transit case (panel (a)) differ 
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Figure 5. QU plots in transit cases (a) and in bypass cases (b): In both cases the source size p* = 10 — 3 and the direction of the relative velocity V is assumed 
to be parallel to the observer's x-axis (a = 0). In the left panel, the ratios of uo to p* are 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2 and 0.1 from the leftmost curve. 
In the right panel, they are obtained as 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9 and 2.0 from the outer curve. 



considerably from those in the bypass case (panel (b)). As long as p* <g 1, however, we can specify the shapes of these curves by one 
parameter, uo/p,, because the polarization degree depends only on the parameter u/p, for small p„. Hence, if we can fit the data of the 
Stokes parameters on an appropriately scaled curve, we can estimate uo/p*. 

Besides the possibility of estimating p», we should note that the QU plot provides information on the direction of the relative microlens 
motion V. Microlensed QU plots are shown in Fig. |6| (panel (a)), in the cases when the microlenses move along lines inclined to the 
observer's a>axis with an angle a — 15° (panel (b)). In these cases, we find that the plots lean at an angle —2a as compared to the ones with 
q = (as shown in Fig.|5J. In Fig.|6| we draw the slanted QU plots for three cases: when the microlens moves from the lower-left corner 
to the right in panel (b) in the bypass case (no = 0.0015), the same movement in the transit case («o = 0.0003) and when it moves from 
the left to the upper-right corner in the bypass case (uo = 0.0012). While, in the first two cases, a point on the QU plot moves clockwise 
around the origin with time, in the last case, the point moves anti-clockwise. Hence, we find that points on the QU plots, in general, move 
anti-clockwise when the vector u x V points to the observer from the lens plane, and vice versa. 

Thus, if we have data on the Stokes parameters, i.e. (Si, Sq, Su) during a microlensing event, we can obtain not only the finite source 
effect p, but also the direction of the relative microlens motion, including the sign of u x V . 



5 SUMMARY AND DISCUSSION 

In this paper we have investigated the single microlens effect on the Stokes parameters in order to obtain some semi-analytical formulae. As 
a result, we observe that the formulae not only reduce the double integrals in the estimations of the parameters but also can be expressed as 
a useful approximation form in the bypass case. In addition, we have shown that by combining the polarimetric data with photometric data, 
we can estimate not only the finite source effect p* but also the direction of the microlens motion, including the sign of it x V. 

Although we have argued on the benefits of observations of polarized stars by the microlensing event, we have not yet enumerated the 
difficulties involved in the polarimetric observation of the microlensed star. In this regard, we shall comment on the following difficulties 
in the remainder of this paper: (1) the lower limit of the detectable polarization degree, (2) the limiting magnitude of the polarimetric 
observation, (3) the observable duration for polarization and (4) the identificati on of the target star . 

The first difficulty is concerned with observational technology. Althoug h lKempetalJil983h succeeded in measuring the polarization 
degree with ~ 0.01 per cent accuracy during an eclipse in a binary system, it is not easy to detect polarized light with such accuracy even 
with the present technology. At present, we can expect the accuracy of the polarization degree to be Pu m ~ 0.05 per cent (Hirata, in private 
communication). This limit defines the upper limit for u as u < 3p* for p„ < 0.1 (see Fig.|3j. Hence, we can define the effective radius R 
in the polarimetric observation of a microlensed star that depends on Py im , e.g. R w 3p* for Py lm = 0.05 per cent. 4 

The second difficulty arises due to the fact that the limiting magnitude of polarimetry is in general 5 — 7 mag brighter than that of 
photometric observation. Most of microlensed stars have magnitudes in the range of 15 — 21 mag in the photometry (/-band). This means 



4 The effe ctive radius, of course, depen ds on coefficients ci , C2 in equations |8j and 1101 . The value estimated here is in the case where c\ = 0.64 and 
c 2 = 0.04 ISchneider & WagoneJ 19871) . 
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Figure 6. QU plots and the corresponding trajectories of the microlenses: When the microlens moves from 1 to 7 on each line in panel (b), the point on the 
QU plot moves from 1 to 7 on the corresponding curve in panel (a). In this figure, it is assumed that p* = 10 -3 and the lines are inclined to the x-axis with 
a = 15°. The time interval between the circled points is 5 X 10 -4 tE- 



that the corresponding magnitudes in the polarimetry of the star are in the range of 20 — 28 mag. Therefore, in order to detect the polarization, 
we need a large telesco pe. In FOCAS of the Subaru Telescope, the limiting magnitude is 26.2 mag in the /-band with a 1200 s exposure 
iKashikawa et all2002l) . Thus, the observation is very difficult even by using FOCAS. Fortunately, while the polarization is detectable, the 
magnification factor of the star is larger than 1 / R. Hence, the polarized light from the microlensed star with even the faintest magnitude may 
be detectable by FOCAS because the magnification factor is larger than 10 when p* < 0.03. 

The third difficulty is associated with the exposure time for the microlensed star. For the observation of faint stars, a long exposure time 
is necessary. In order to extract the finite source effect, we need to observe the lensed star a number of times. Further, the total exposure 
time must be less than the observable duration for polarization, which is equal to the crossing time At c of the microlens over the effective 
diameter (= 272), i.e. 2Rt^ m 6p,te- We thus obtain the lower limit of p» in the polarimetric observation of the microlensed star as follows: 



4.6 x 10" 



Not 



100 



1200 s 



50 days 



< P*, 



(27) 



where N b s and t cxp denote the frequency of the observations and the exposure time, respectively. 

Among the four difficulties listed above, the last one may be the most challenging. It arises from a reason specific to the microlensing 
event. Generally, microlensed stars are observed in very crowded regions such as our Galactic bulge, LMC and SMC. Hence, there are many 
stars around the target star within the same field. As 50 — 100 stars can be simultaneously observed using FOCAS (Kawabata, in private 
communication), it may be possible to detect the polarization. 

If we can overcome the above difficulties, we can estimate the fi nite source effect from the polarimetric observation and decouple the 
degeneracy in the microlens parameters by using the relation given bv lvan Belle! Il999l) . Even after satellite projects succeed in measuring 
the angular Einstein ring radius 0e, the polarimetric observation of the microlensed star will provide us eful information such as on the limb- 
darkening model, on the relation between angular diameter and colour-magnitude of stars Ivan Belle!) and so on. We thus believe that it is 
important and w orthwhile to fo cus attention on polarimetry in microlensing events after receiving microlens alerts such as from the Early 
Warning System ludalskikoojh . 5 



http : / / www . astrouw .eclu.pl/~ogle/ ogle3/ews/ews . html 
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APPENDIX A: ELLIPTICAL INTEGRALS 

The elliptical functions of the first, second and third kind (F, E, If) are given as follows: 




(Al) 



(A2) 



- 1 







Jo (l-nsm 2 9)y / l-k 2 sm 2 6 
Let us define the integration as follows: 
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Jo 



(1 - nt 2 )^/(l -t 2 )(l - k 2 t 2 )' 



dt 



(A3) 




(A4) 
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The objective in this appendix is to express A; as a combination of F, E and II. First, we consider the derivative of \tyj f(t)]/(l — nt 2 ) 1 1 
with respect to t, where f(t) = (1 - t 2 )(l - k 2 t 2 ): 



d_ 

dt 



(1 - nt 2 ) 1 



2(1 - l)/(t) + (1/2) (1 - nt 2 ) {tf(t) + 2(3 - 2J)/(t)} 



(A5) 



We restrict ourselves to I ^ 2. From the definition of A;, the first three functions, Ai , A , A_i, are expressed by the following three complete 
elliptical functions, 

Ai(t;/c,n) = n(sin _1 t;fc,n), A (t; k, n) = F(sin _1 t; k) and A_i(t; fc, n) = (1 - ■^■)_F(sin" 1 1; k) + ^E(sin~ 1 1; k), (A6) 
respectively. We denote the numerator of equation <A5> as g(t) 

g(t) = 2(1 - l)/(t) + -(1 - nt 2 ) [tf'(t) - 2(21 - 3)f(t)] =A + B(1- nt 2 ) + C(l - nt 2 ) 2 + D(l - nt 2 ) 3 , (A7) 
where 



1 



A = 2(1 - 1) 1 - - 1 , B = -(21 - 3) 1 - 2 



1 + k 2 



+ 3^r , C = 2(Z-2) 3— 



1 + k' 



andD = -(21 - 5) — . 

n 2 



These coefficients enable us rewrite equation <A5> as follows 

t 
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dt 



(1-nt 2 ) 1 



dt dt dt dt 
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+2(1-2) 3 
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dA ; 
n j dt 

1 + fc 2 \ dA ; 



dt 



- (21 - 3) 1 - 2 



(2/ - 5) 



,1 + fc 2 
n 

W-3 



fc 2 \ dA,_! 



dt 



(A8) 



We integrate the above equation from t = to 1 in order to obtain the following recursion: 

2(l-l)(l--) ( 1-— \ A,(l;*,n) = (21-3) f 1 - 2^+^+3^) Anfl;^) 
V n J V n / y n n J 



-2(1-2) 3 



1 + k' 



A ! _ 2 (l; fc, n) + (21 - 5)- ? A i _ 3 (l; fc, n). 



For example, in the cases of I = 2, 3, we find 



2 1 



1 



A 2 (l;fc,n) 



1-2 



1 + k 3k 2 



n(7r/2;fc,n) + 



F(n/2;k)- -E(n/2;k), 
n 



(A9) 



(A10) 



and 



1 + k 2 



:ik- 



1 + k' 



4(l--j ^-_jA3(l; fe; n)=3^1-2-- + -jA 2 (l; fe) n)-2^-- 

We substitute n = 4pu/(u + p) 2 ,k 2 = 4n/[(u — p) 2 + 4] into equation <A9t to obtain a recursive formula as follows 

'l-2\ (p 2 - u 2 ) 2 - 8(p 2 + u 2 ) 



n(7r/2; k, n) + -%F(n/2; fc).(All) 

nr 



Ai(l;k,n) 
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(u- p) 2 



I- 1 



(p-uY 



-A { _ 2 (1; k,n) 



21 -5\ (p + u) s 



/2/ - 5\ 



A;_ 3 (l; fc, n). 



For example, in the cases of I = 2, 3, we obtain 

A 2 (l;fc,n) = " § + p3 ~ 2 n(V2;fe,n)-l F(7r/2;fc) + 

(u- p) 2 2 \u- pj 

am, ^ M 4 + 4u 2 / 9 2 -2(u 2 + p 2 )+p 4 + 6 TT , .„ , , 
A 3 (l;fc,n) = — II(7r/2;fc, n) 



(u + p) 2 {(u-p) 2 +4) 
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4(u - p) 6 
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E(n/2;k). 
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